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Propagation Speed of Gravity and the Relativistic Time Delay 



We calculate the delay in the propagation of a light signal past a massive body 
that moves with speed v , under the assumption that the speed of propagation 
of the gravitational interaction c g differs from that of light. Using the post- 
Newtonian approximation, we consider an expansion in powers oiv/c beyond the 
leading "Shapiro" time delay effect, while working to first order only in Gm/c 2 , 
and show that the altered propagation speed of the gravitational signal has no 
effect whatsoever on the time delay to first order inv/c beyond the leading term, 
although it will have an effect to second and higher order. We show that the only 
other possible effects of an altered speed c g at this order arise from a modification 
of the parametrized post-Newtonian (PPN) coefficient a\ of the metric from the 
value zero predicted by general relativity. Current solar-system measurements 
already provide tight bounds on such a modification. We conclude that recent 
measurements of the propagation of radio signals past Jupiter are sensitive to cti, 
but are not directly sensitive to the speed of propagation of gravity. 

Subject headings: gravitation: experimental 



The time delay of light, also known as the Shapiro effect, is one of the observational 
cornerstones of general relativity. For a light ray passing a body of mass m a at location x a , 
and received by an observer at x r , the excess delay is given by (modulo a constant) 
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where k is a unit vector in the direction of the incoming light ray, x TO = x r — x a , and 
r ra = |x m |. If the ray passes by the body a minimum distance d that is small compared to 
r ra , then 

2Gm a (2r ra \ 

The one-way delay amounts to 70 microseconds for a ray that grazes the Sun and 10 nanosec- 
onds for a ray that grazes Jupiter. In addition to tests of the time delay by the sun (Reasen- 
berg et al. 1979), the time delay has also been studied in a number of binary pulsar systems 
(see for example, Stairs et al. (1998)). 

Advances in precision timekeeping and in VLBI have motivated the study of the first 
relativistic corrections, of order v/c, to the basic Shapiro formula. These corrections would 
be relevant for signals passing near moving bodies, such as bodies in binary pulsar systems, 
or planets such as Jupiter. In a series of papers, Kopeikin and colleagues have analyzed these 
correction effects in detail (see Kopeikin and Schafer (1999) and references therein). They 
find that, to first order in v a /c, where v a is the velocity of body a, the corrected expression 
takes the form 

2Gm a f k • v a \ i 
A = — I 1 1 ln[r ra - K • x ra J 

f (3) 

where 

K = k - -k x (v a x k) , (4) 

c 

where all quantities are to be evaluated at the same moment of time, say, the time of reception 
of the ray. 

For a light ray passing close to Jupiter, these corrections are small, of order picoseconds 
or less, but may be detectable with the latest VLBI techniques (Kopeikin and Fomalont 
2002). 

Recently, however Kopeikin has argued that, because of the motion of a body such as 
Jupiter during the passage of the light ray, the fact that the gravitational interaction is 
not instantaneous should have an effect on the time delay (Kopeikin 2001). To study this 
phenomenon, he allowed the speed of gravity to be c g , where c g ^ c; this would modify the 
retardation of the gravitational fields used to calculate the time delay. He argued that Eq. 
(3) would be modified by replacing c by c g in the prefactor 1 — k • v a / c in Eq. (3), and in the 
velocity-dependent term in Eq. (4) (Kopeikin 2002). Kopeikin and Fomalont (2002) then 
proposed a test of the propagation speed of gravity via VLBI measurements in the fall of 
2002. 
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This result goes against intuition. First, while the retardation of gravity surely has 
effects, one does not expect to find effects odd in powers of v/c until the leading effects 
of gravitational radiation reaction occur, far beyond the order of the -u-dependent effects 
claimed by Kopeikin. Second, the time delay is controlled dominantly by the Newtonian 
gravitational potential U; if this is a retarded potential with speed c g , it would be given by 

^.^■/ ^-'pJ-^ dV, (5) 

where p is the density of matter. Expanding the retardation in the near zone (|x — x'| <C 
a characteristic wavelength of gravitational waves), where the time delay predominantly 
occurs, we obtain 

J |x - x'| c g dtJ I K ' J 

+ 2^^/ p(t ' X ' )|x - X,|rfV + ---' (6) 



The second term would have been of order v/c g relative to the leading term, and might 
have contributed to the time delay, but it vanishes because of the conservation of mass. The 
effect of c g first appears in the third term, which is of order (v/c g ) 2 relative to the first. This, 
rather than first order in v/c g , is what intuition dictates would be the first level at which 
retardation-dependent effects occur. This argument was invoked implicitly in the paper of 
Asada (2002). 

Another argument against such v/c g dependent effects comes from electrodynamics. It 
is a well-known fact that the electric field of a uniformly moving charge points not toward 
the retarded position of the charge, but toward the present position (see, for example, §11.10 
and Exercise 11.17 of Jackson (1999)). In gravity, because acceleration involves Gm/r ~ v 2 , 
working to first order in v/c is equivalent to assuming a uniformly moving body, and thus a 
similar conclusion holds: to order v/c, retardation should have no observable effects. 1 

These considerations have motivated us to check Kopeikin's claim with a careful calcu- 
lation of the time delay to the same order. We have assumed that, to the post-Newtonian 
order required for the calculation, the metric can be described by a matter gravitational 
potential U and a current gravitational potential V 1 , both of which are assumed to be re- 
tarded potentials using a speed c g , as in Eq. (5). The terms in the metric also are multiplied 
by suitable parametrized post-Newtonian (PPN) parameters to reflect the possibility that, 



1 A similar conclusion was reached by Carlip (2000) in the context of binary motion, in response to a 
proposed bound on the speed of gravity by Van Flandern (1998). 
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whatever theory one is considering with a different speed of gravity, it will necessarily not 
be general relativity, and thus may have different PPN parameter values (see Will (1993) for 
a review of the PPN framework). In such a theory, those PPN parameters may well depend 
on c g , but in the absence of a specific theory, we leave them general. 

We then calculate the time delay to order (Gm a /c 3 )(v a /c), for a signal emitted from 
a distant source and received on Earth, passing through the gravitational field of the solar 
system. We carefully treat the propagation of the light ray both through the far zone, where 
the propagation time is long compared to the dynamical time of the system, and through 
the near zone, where the propagation time is short compared to the dynamical time, and 
where positions and velocities can be expanded about values at some common time, such as 
the time of reception on Earth. Our result closely parallels Eq. (3): 



with K given by Eq. (4). 

The first important observation is that c g appears nowhere in this expression. Although 
it appears in various intermediate expressions, it cancels in the end, to the order of ap- 
proximation shown. In addition to the overall factor of (1 + 7)/2, which is standard, the 
only change comes in the coefficient (1 + () in front of the velocity-dependent term in the 
prefactor. Here ( is a function of the PPN parameters ct\ and 7, given by 



In general relativity, 7 = 1, a± — 0, and we recover Kopeikin's result, Eq. (3). 

Secondly, the velocity dependent corrections in Eq. (7) have a simple physical inter- 
pretation. Since we are working to first order in v/c, only Jupiter's rectilinear motion is 
relevant; its acceleration generates effects higher order in Gm/c 2 . The vector K is simply 
the spatial unit vector of the photon's direction as seen in the rest frame of Jupiter, the 
correction terms in Eq. (4) are nothing but the effect of aberration. They depend only on c, 
which is the velocity relevant for Lorentz transformations and the propagation of light, not 
on c g . Indeed, from the viewpoint of Jupiter's rest frame, c g is completely irrelevant, since the 
gravitational field is static (again, ignoring Jupiter's acceleration). The velocity- dependent 
prefactor is a combination of a gravitomagnetic effect (the fact that the gravitational field is 
not a scalar quantity, but contains both a vector and tensorial part) and a simple Doppler 
effect in transforming from Jupiter's rest frame to the barycentric frame. 
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We thus conclude that the v/c corrections to the Shapiro time delay are normal 1.5PN 
corrections that occur when there are moving bodies, but that they have nothing to do 
with the speed of propagation of gravity, insofar as it affects the retardation of gravitational 
interactions. Furthermore, as a potential test of alternative gravitational theories, measuring 
these v/c terms is not promising, because a variety of solar system measurements already 
constrain a.\ and 7 to such a degree that \(\ < 4 x 1CT 3 under relatively weak assumptions 
or \(\ < 5 x 1CT 5 under assumptions that invoke so-called "preferred-frame" tests of the 
parameter cti (see Will (2001) for the latest bounds on the PPN parameters). In fact, the 
VLBI measurements are sensitive mainly to the velocity dependence in the logarithmic term, 
not to the prefactor. Therefore, measurements of the propagation of radio waves past Jupiter 
do not directly constrain the propagation speed of the gravitational interaction. 

The remainder of this paper provides details. In Sec. 2 we describe the assumptions 
and basic equations that go into our calculation. Section 3 carries out the integrations to 
find the time delay. In Sec. 4 we give concluding remarks. 



We assume that, whatever theory of gravity is in force, it is a metric theory, that is it 
has a spacetime metric g^ v that governs the interactions and motions of all non-gravitational 
fields and all "test" particles. In particular, we will assume that light rays move on null 
geodesies of the spacetime metric. We will assume that atomic clocks measure proper time 
as given by the invariant interval of the spacetime metric. By choosing the units determined 
by such physical measurements appropriately, we may make the speed of light as measured 
by any freely falling observer unity. Henceforth we will use units in which G = c = 1. 

We then assume that the theory of gravity has a post-Newtonian limit that can be 
written in the following form: 



2. Assumptions and Basic Equations 



9oo 



l + 2f/ + 0(e 2 ), 

2(1 + 1 + ai /A)Vi + 0(e 5/2 ), 




(9) 



where U and Vi are retarded gravitational potentials given by 
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where p and v l are the density and velocity of matter, and where we have inserted c g in the 
retardation expression to reflect the possible difference between the speed of propagation of 
gravity and that of light (unity). The parameter e is a bookkeeping parameter inserted to 
keep track of "orders of smallness" of various quantities, with v ~ e 1//2 , U ~ e, ~ e 3 / 2 , and 
so on. 

To the orders of approximation shown, the metric in Eq. (9) is general enough to 
encompass a broad class of metric theories of gravity with a propagation speed of gravity 
different from unity, with the exception of "massive graviton" theories (Visser 1998; Babak 
and Grishchuk 2002), whose gravitational potentials also include Yukawa-like modifications. 
In the special case of 7 = 1, ct\ — and c g — 1, the metric corresponds precisely to that of 
linearized general relativity in harmonic gauge, to the order considered. There are two free 
PPN coefficients 7 and cti, the coefficient of U in g 00 has been chosen to be unity as usual 
by our choice of the measured value of the gravitational constant. We ignore all "preferred 
frame" terms that could appear in the PPN metric (see Will (1993) for discussion). 

Notice that the terms ignored in each metric component are all 0(e) larger than the 
terms kept; no correction terms a "half-order" higher are present, such as 0(e 3 / 2 ) terms in 
goo or gij, or 0(e 2 ) terms in g 0i . This has certainly been true for all those theories analysed by 
the present author or by others to date, including several, such as the Rosen bimetric theory, 
with c g 7^ 1 (Lee et al. 1976). Such correction terms generically arise from non-linear terms 
in the field equations, which are minimally of order m/r ~ e higher than the leading term, 
or from 0(v 2 ) ~ 0(e) corrections in the definition of the stress-energy tensor. However, the 
retardation effects in the potentials could in principle lead to such "half-order" terms in a 
slow-motion expansion. It is precisely these terms whose effects we wish to explore. 

Before proceeding, however, we emphasize that we have assumed a specific theoretical 
framework for the analysis. Relatively few theories with c g 7^ c have been analysed within 
the PPN framework (the Rosen theory and a few others). Furthermore, it is entirely possible 
that c g ^ c theories may be more fully or consistently realizable in the context of non-metric 
theories, for which the PPN framework is not applicable. In the absence of a concrete 
example, we will restrict attention to the PPN approach described above. 

In a metric theory of gravity, the gravitationally coupled Maxwell equations require that 
light rays propagate on null geodesies of the metric (9), given by the equations 



where = dx^/dt (N° = 1) along the light ray. 

Calculation of the time delay follows the method of Will (1993): write the trajectory of 



d 2 x i /dt 2 + N^N u {T 1 ^ - 
(ds/dt) 2 = g 00 + 2^ + gijrfN* 



0, 
0, 



(11) 
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the light ray in the form 

x l (t)=x l e + k l (t-t e )+x l p (t), (12) 

where e denotes the event of emission of the signal, k l is a spatial, Cartesian unit vector, 
and x p (t) denotes the perturbation in the path of the ray, which is clearly of order e. Then, 
substituting N l = k l + x % and = i]^ + h^ u , where rj^ is the Minkowski metric, into the 
second of Eqs. (11), we obtain k • x p = — h^ u k^k u /2 + 0(e 2 ), where fc M = (1, k l ) and overdot 
denotes d/dt. But from Eq. (12), we see that |x — x e | = t — t e + k • x p + 0[\x p \ 2 /(t — t e )], 
consequently, 

t-t e = |x - x e | + A(t,t e ) , 

i r* 

A(Me) = -FF / V[t,x(t)]rft + 0(e 2 )|x-x e |. (13) 

Equation (13) is completely equivalent to Eqs. (20) and (44) of Kopeikin and Schafer (1999) 
and to Eq. (18) of Kopeikin (2002). To the required order, the function x(t) that is used to 
evaluate in the integral is the unperturbed path of the ray, given by x(i) = x e + k(t — t e ). 

We now specialize the metric to the physical situation of a set of bodies of mass m a , 
and velocity v a (the solar system). We work in a coordinate system whose origin is at the 
barycenter of the system, so that J2 a m a v a = 0. We parametrize the trajectory of each 
body by x a (w), where u is a parameter that is proportional to coordinate time (with unit 
proportionality constant), such that x a (w = t ) is the location of the body at coordinate 
time t = to, x a (-u = 0) is the location of the body at coordinate time t — 0, and so on. 
The ordinary velocity is given by v a (u) = d?t a (u)/du = dx. a (u)/dt. (We introduce this 
parametrization merely to avoid confusion with other "times" to be discussed.) Treating 
each body as effectively a "point" mass, with p(t',x') = X] a m «^ 3 [ x ' — x a("" = t')] we find 
that the potentials U and Vi take the Lienard-Wiechert form, 



= En- 



m a 



x - Xo(s )| - v a (s a ) • (x - x a (s a ))/c 9 ' 



W' x ) = El ? m m f[ Sa ) nw-> (14) 

^ |x-x a (s a )| - v a (s ; • (x-x a (s a ))/c 9 

where s a is "retarded" time given by the implicit equation 

s a = t-\x-x a (s a )\/c g . (15) 

It is important to understand precisely the meaning of the expressions x a (s a ) and v a (s a ). 
They do not mean that x a and v a have become functions of s a ; they are functions only of 
the parameter u. Rather the expression x a (s a ) means x a (w = s a ), or "x a (u) evaluated at 
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u = s a where s a is given by evaluating Eq. (15) for a chosen field point (t,x)". This will be 
important, for example, when we carry out Taylor expansions of x a (s a ) and v a (s a ). These 
issues are addressed in detail in an Appendix. 

Substituting Eqs. (14) into Eqs. (9) and thence into (13), we can write the time delay 
in the form 

m, t\ r,x iV /"' [l-(2 + Qk-v,,(s„)]<« 

A(t„, t.) = (1 +7) E m. I | x _„.(,.), _ v .(,.) .(,_„.(,.))/,, • («) 

where t r denotes the time of reception of the ray and where £ = ot\/{2 + 27). 



3. Calculation of the time delay to 1.5PN order 

Equation (16) is a sum of contributions for each body in the system. Since we are 
working to linear order in to, we can evaluate the integral for each body separately, then 
multiply each by m a and sum over the bodies. We rewrite the unperturbed trajectory of the 
light ray in the more convenient form 

x(*)=k<7 + £ , (17) 

where a is the time t modulo a constant, and £ a is a vector from the barycenter to the light 
ray at the moment of its closest approach to body a; the parameter a is chosen so that a = 
at this point. Because we are working to first order in m a , we can ignore all effects related 
to the deflection of this ray. 

We adopt the simplified notation 

T a (<r, s a ) = x(o-) - x a (s a ) , 
r a {<r,s a ) = |x(o-) - x a (s a )| , 

s a = a-r a (a,s a )/c g , (18) 

where we have dropped the irrelevant arbitrary constant that initializes time. The integral 
to be evaluated is 

* M =r ,"-f?, T ;"-r m 

Ja e r a (a, s ) - -v a (s ) • r a (a, s a ) 

We now divide the integral into two contributions ^i 1 ^ and , where ^i 1 ^ = JJ 1 . . .da 
and ^i 2 -* = . . .da where the time parameter o\ is chosen to be large and negative (recall 
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we chose a = at closest approach), and to satisfy the following inequalities |x a | <C |o"i| <C 
(l/f a )|x a |. Since v a <C 1, an appropriate value of <j\ clearly exists. Then in evaluating ^i 1 " 1 , 
we can expand the integrand in powers of 1/cr, since \a\ ^> |x a | everywhere, and integrate 

(2) 

by parts. In evaluating ty a , since the integration is always within the near zone, we can 
expand the variables of the bodies about their values at the coordinate time corresponding 
to closest approach in a slow- motion expansion. We will work only to first order in v a . Each 
integral will depend on the arbitrarily chosen value of U\. We then match the integrals by 
expanding the relevant terms in in a slow-motion expansion, and the relevant terms in 
$ a 2 ^ in powers of \jo\. We will see that all dependence on o\ cancels to the order considered, 
as it must. 

We first evaluate \I/i 2 ' '. We re-initialize the parameter u for body a so that it is zero at 
the time of closest approach and expand the variables of body a about that point, so that 
Xo(ii) = x a (0) + v a (0)w + 0(eu). When evaluated at that value of u equal to retarded time 
s a , we get an analogous result x a (w = s a ) = x a (0) + v a (0)s a + 0(es a ). Also, v a (u = s a ) = 
v o (0) + O^Va). Note that x o (0) and v a (0) are evaluated at u — 0. By defining the vector 
d a (s a ) = £ a — x a (s a ), which points from the body at retarded time to the point of closest 
approach, we can write r a (a, s a ) = ker + d a (s a ). We note that, to the necessary order, the 
denominator in Eq. (19) can be written r a (a, s a ) — v a (s a ) - r a (a, s a )/c g = |D| +0(er a ), where 

D = r (<r, s a ) -r a (a, s a )v a (s a )/c g 
= ka + d a (s a ) - v (s )(<r- s a ) 
= ka + d o (0) - v o (0)a + [d a (s a ) - d a (0) + v a (0)s a ] 

+0(er a ) . (20) 

The term in square brackets in Eq. (20) vanishes to the necessary order. Notice that c g has 
cancelled out. 

Defining the vector z(a) = ka + d a (0), with k • d o (0) = (k and d a are orthogonal at 
the moment of closest approach to body a), we can write D = z(er) — v a (0)<7, and thus 

r a {cr,s a ) -v a (s a ) ■r a (<j,s a )/c g = |z (a) - v a (0)a\ + 0{er a ) 

= [a 2 (l - 2k • v a (0)) + 2<j(k - v o (0)) • d a (0) 

+4(0) 2 ] 1/2 + O(er a ). (21) 

Inverting this expression and integrating with respect to a yields the result 

*i 2) = [1 " (1 + C)k • v a (0)] In { (z{&) + k ■ z(a)) (1 - z(a) ■ v o (0))}^ + 0(e) . (22) 

At the reception point, o>, we define z(a r ) = ka r + £ — x a (0) = x r (o>) — x a (0) = St ra and 
notice the fact that <i a (0) 2 = f 2 ra — (k • x ra ) 2 . At <7i, we expand the integral in powers of 
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d a /Wi \ «C 1, using the expansions (to the required order) 

z(<7i) = -<ti - d a (0) 2 /2<n + 4(0) 4 /8^ 3 + O(d a (0)Vr 5 ) , 



= -(4(0) 2 /2<Ti)(l - 4(0) 2 /4<t 2 ) + 0(4(0) V 5 ) , 

-k-d a (0)M + kK(0) 2 /2 ( T 2 ) + O(rf a (0)Vr 3 ). (23) 

Substituting these expansions into Eq. (22), and keeping terms through 0(crf 2 ), we obtain 



z(a 1 ) + k • z(<7i) 

Z((7i) 



+ 



l-(l + C)k-v o (0) 
l-(l + C)k-v o (0) 



i ^ra k • X ra 

IU| 2 

4(0) 2 



X ra • V o (0) 



-k.V o (0) 



In I ail + 



d.(0).v.(0) + (K^ + 0W())V) 



0-1 



Next, we evaluate Here, because a is integrated over a range that could include 

many dynamical periods of the masses as the signal propagates from the distant star toward 
the near zone, we must leave all variables such as x a (s a ) unexpanded. Instead, we expand 
all expressions in advance in powers of I /a. Then 



(24) 



r a (cr, s ) = k(T + d a (s a ), 
r a {cr,Sa) : 



, , f v 4(Sa) 2 . k- d a (s a )rf ± (s a ) 2 4 _ 3 

-cr - k • d a (s a ) 7T — + — + 0{d a {s a ) a 6 ) , 



2a 



1 + - (T + 



k-d a (s a ) c/±(s a ) 2 k • d a (s a )d±(s 



+ 



2Cg<7 



2c g a 2 



+0(d a (s a ) 4 a- 3 ) , 



(25) 



where d±(s a ) = d a (s a ) 2 - (k • d a (s a )) 2 . Calculating [r a (<7,s a ) - v(s a ) ■ r a (cr, s a )/c g ] _1 , and 
expanding it in powers of a^ 1 , and then expanding to first order in e 1//2 , we obtain the 
integral 



- (k ■ d a (s a ) + d "( g ") " Va ( s ") _ 2r]l k-d a {s a )k-v a (s a ) 

2 \ r- n 



a 

-(- 

a 3 12 



K( So ) 2 -3(k.d ( Sa )) 2 ]-— (So;r '" 2 5 - , '-^ 21 



-|^ a ( Sa ) 2 -V(k-d a ( Sa )) 2 ] 



-2 



k • d a (s a )k • v a (s a ; 



+ 



^ (3rf a (s a ) 2 k • d a (s a ) - 5(k • d a (s a )) 3 ) } 



+0(e,rf a ( Sa )V" 4 ), 



(26) 
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where 77 = 1 + (2 + ()c g /c, 7/ = (1 +r))/2 and 77" = (5 + 3r])/8. We then integrate the various 
terms by parts, using the formula f(da/a n+1 )f(s a ) = —f(s a )/na n + f(da/na n )df /da. Each 
cr-derivative df /da raises the order of that term by e 1 / 2 ; we stop integrating by parts when 
the resulting term is 0(e) relative to the leading term. Because a in the integral parametrizes 
the light path, and because we are not Taylor expanding the orbital variables, we must take 
into account that s a depends on a through Eq. (25), so that 



d 

da 



d a (s a ) = -Vo(So) 



ds a 
da 



-y a {s a 



1 + 



1 



d±(Sg) 21 

2c 9 a* 



+ 0(e,d a (s a ) 3 a- 3 ) 



(27) 



After integration, we discard terms that depend on negative powers of the emission time a e ; 
for light from a distant quasar, these are clearly negligible. The result is 



l-(l + C)k-v a ( Sl ) 



1 

In \a 1 \ + - 

c 



l-(l + C)k-v a (s e 



In \a P . 



--k • d a ( Sl ) - — K( Sl ) 2 - 3(k • d a ( Sl )) 2 ] 



1 



6af 
1 



4a( 

[3rf a ( Sl ) 2 k-d a ( Sl )-5(k-d a ( Sl )) ; 



2a! 
1 



Aa\ 



1 - - d (si) • V a ( Sl ) + 1 + 2C + - k • da(si)k • V (si) 



2 (1 - 1^ k • d a ( Sl )d a ( Sl ) • v (si) - (1 + ()d a ( Sl ) 2 k • v (si) 



1 



+ 1 + 3C + 2- k-v a (si)(k-d ( Sl ))' 



+ 0(e,d a (s a ) A a^). 



(28) 



Since <7i is within the near zone, we can use the expansions 
d a ( Sl ) = d o (0) - v a (0)si + 0(ed a ) 



d a (0)-v a (0) 
v a (0) + O(e), 



1 + 



°i + 



4(0) 
2c„ai 



21 



+ O(e,rf a (0)V 3 ) 



(29) 



where, again, (0) means w = 0. Recalling that k • d a (0) = 0, and substituting these expres- 
sions into Eq. (28), we obtain 

1 



1 - (1 + C)k • v a (s e ) In \a e \ + k • v a (0) 1 + 



l-(l + C)k-v a (0) 



In |<7i| — 



dq(0f 

Aa\ 



d o (0)-v a (0) (l-C\ 4(0) 2 k-v a (0) 



+ O(e,rf a (0)V')- 



(30) 
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Combining the expressions (24) and (30) for and ^i^, we see that all terms involving 
o\ cancel, as expected. Our final expression for \l/ a is 



+ 



l-(l + C)k-v a (0) 
l-(l + C)k-v a (s c ) 



In 



In \ a P 



2 

0(e). 



x m -v a (0) k-v a (0) 



(31) 



Notice that the dependence on c g in the term k • v a (0)/c s in Eq. (31) is illusory: to 
obtain the measured effect, we must multiply \l/ a for body a by its mass m a , and sum over 
all the bodies in the system. Although the event corresponding to a = (closest approach) 
is different for each body, the effect of these differences on v a is of order of the acceleration 
of each body, which is 0(e), hence we can evaluate v a (0) at the same coordinate time for 
each body, say, the time of reception of the ray. But, because 'Yl ia Tn a v a = in barycentric 
coordinates, the net effect of this c 3 -dependent term vanishes. Similarly, for a given light 
ray, the time a e is a constant, independent of the body. The difference between the various 
values of s e for different bodies has an effect only at order e, so again, summing over all the 
bodies in the system causes the velocity correction to the In \ a e \ term to vanish. The In \ ca e \ 
term then is an irrelevant constant. By the same argument, the factor of 2 in the logarithm 
produces only an irrelevant constant. The final formula for the time delay, correct to 1.5PN 
order is 



A(t r , t e ) = -(l + 7 )5> a {[l-(l + C)k- v o (0) 



In (f r 



k ■ x ra ) 



'q(O) | 

n ' 



• (32) 



The speed of propagation of gravity, which initially appeared in the retarded potentials, is 
nowhere to be found in this final expression. The only consequence of using an alternative 
theory of gravity is in the PPN parameters 7 and £ = cti/(2 + 27). 

The variable x m = x(cr r ) — x a (0) can be converted so that all times are referred to the 
moment of reception of the ray, using the expressions 



^-ra X ra -|- V a CT. 

= X 



v o (0) = v a + 0(e), 



(33) 



where x m = x(o>) — x a (cr r ) and v a = v a (ay). Substituting into Eq. (32) and rearranging, 
we find, to the necessary order, 



A(t r , t e ) = -(1 + 7) J2 m *{ C)k • v tt ] In (r ra - K • x ra )} , 



(34) 
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where 

K = k — k x (v x k) . (35) 

In general relativity (7 = 1, £ = 0) this agrees completely with results derived by Kopeikin 
and Schafer (1999), to 1.5PN order. 

4. Discussion 

We have shown that the speed of propagation of gravity has no direct influence on the 
time delay of light to 1.5PN order. The only effect comes from any modification of the PPN 
parameters that might arise in a theory with a different propagation speed. This contradicts 
claims made by Kopeikin (2001, 2002). 

In addition, existing solar-system experiments already place sufficiently strong bounds 
on the PPN parameters 7 and a.\ that even the potential difference from GR in the 1.5PN 
correction represented by the parameter ( in Eq. (34) must be small. The parameter 
7 is known to be unity to a part in 1000 from standard time delay and light deflection 
measurements, while \ai\ < 2 x 10~ 4 from analyses of Lunar laser ranging (Miiller et al. 
1996) and binary pulsar data (Bell et al. 1996). The latter bounds assume that the solar 
system and the relevant binary pulsar are moving with respect to the cosmic background 
radiation with known velocities of order 300 km/s, and that, in any theory of gravity with 
«i 7^ 0, the rest frame of that radiation coincides with a cosmological preferred frame. Such 
a frame is to be expected in a theory of gravity in which c g 7^ c for the following reason: in 
a theory with either interacting dynamical fields or with non-dynamical "prior-geometrical" 
fields, the speed of gravity is presumably a function of some cosmological values of the fields 
in the theory. But because c g 7^ c, its value depends on the velocity of the frame in which 
it is measured. There must therefore exist some frame in which c g takes its value directly 
from the cosmologically induced values of the underlying fields - the only logical frame is 
the mean rest frame of the universe, and hence the rest frame of the CBR. Then, if ot\ 7^ 0, 
there will be measurable effects in any system that moves relative to this preferred frame. 
With these strong bounds, |C| < 5 x 10~ 5 . 

Nevertheless a weaker bound on a± can be inferred from experiments without appealing 
to preferred-frame effects. We assume only that the theory is a Lagrangian-based theory 
of gravity, so that it possesses suitable conservation laws for total momentum and energy. 
Then, Lunar laser ranging tests of the Nordvedt effect bound the combination |4/3 — 7 — 
3 — 10£/3 — o;i — 2a 2 /3| < 10~ 3 . The perihelion advance of Mercury combined with the 
existing bounds on 7 yield \f3 — 1| < 3 x 10 -3 ; Earth tide measurements force the bound 
|£| < 10~ 3 . One is then left with the relatively generous bound \a± — 2a 2 /3| < 1.6 x 10~ 2 . All 
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bounds on a 2 come from considerations of preferred-frame effects. Thus without appealing 
to preferred-frame tests, one would have to have an extraordinary cancellation between a\ 
and a 2 to conform to Nordtvedt effect measurements, while still making ( large enough to 
make a measurable difference in Eq. (34). 

We are grateful to Steve Carlip, Sergei Kopeikin, Gerhard Schafer, Kenneth Nordtvedt, 
and Hideki Asada for useful discussions. This work is supported in part by the National 
Science Foundation, under grant no. PHY00-96522. 

A. Meaning of retarded time 

Throughout these calculations, quantities such as x a (s a ) and v a (s a ) appear, where s a 
is retarded time evaluated using Eq. (15). In the near zone, we use the expansion x a (s a ) = 
x a(0) + v a (0)s a + . . .. One is tempted to state that the (0) that appears means something 
other than coordinate time equals zero. But this would not be correct. The trajectory 
of the body x a (tt) is a function only of the parameter u, which is chosen to be directly 
proportional to coordinate time. By contrast, s a is not an independent variable, instead it 
is a two-point dependent variable, determined by both the field point (t, x) and the body's 
trajectory. Equation (15) serves to determine (albeit implicitly or iteratively) that value of 
the parameter u = s a at which to evaluate x a (w). The equation quoted above is then simply 
a special case of the Taylor expanded equation x a (w) = x a (0) + v a (0)w + . . ., evaluated at 
u = s a . Figure 1 illustrates this explicitly. On the t', x' plane, assuming flat spacetime, the 
world line of a body moving with fixed velocity v is shown, together with a field point at 
(t, x). The point on the trajectory that corresponds to retarded time from the field point 
is given by the intersection V between the particle trajectory and the past null ray from 
the field point (we assume for this illustration that the speed of the null ray is unity). This 
occurs at a time \x — x a (s a )\ earlier than the field time t, therefore it occurs at a time 
t' — t — \x — x a (s a )\ = s a , which thus corresponds to the value u = s a along the trajectory. 
By considering the projection of this event onto the rr'-axis, together with the location and 
velocity of the particle at t' = 0, it is clear that x(u = s a ) = x(0) + v(0)s a . 
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